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Abstract 


A measurement of the ratio of the branching fractions of the Bg meson to J/i/?fo(980) 
and to ]/tp(p{1020) is presented. The ]/xp, fo(980), and (p{1020) are observed through 
their decays to n^rc^, and K+K^, respectively. The fo and the (p are identified 

by requiring — 974MeV| < 50MeV and |Mk+k- “ 1020MeV| < lOMeV. The 

analysis is based on a data sample of pp collisions at a centre-of-mass energy of 7 TeV, 
collected by the CMS experiment at the LHC, corresponding to an integrated lumi¬ 
nosity of 5.3 fb^^. The measured ratio is g(^2K+K-)^ ~ 0.140 ± 0.008 (stat) ± 

0.023 (syst), where the first uncertainty is statistical and the second is systematic. 


Published in Physics Letters B as doi : 10.1016/ j. physletb . 2016. 02. 04 7. 


© 2016 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license 


*See Appendix|^for the list of collaboration members 









1 


1 Introduction 

Since the observation of the decay Bg —> J/t/^ with J/t/? —)■ , and the mass spec¬ 

trum indicating a large fo(980) component Hi, this channel has been regarded with great inter¬ 
est in heavy-flavor physics. More detailed studies of the 7r+ mass spectrum have shown the 
71+TT^ system to be almost entirely CP odd |l2l|3l. This opens up the possibility of directly mea¬ 
suring the lifetime of the CP-odd part of the meson jUISl. In addition, the Bg —> J/t/^ 71+ 
decay has been used for the measurement of the CP-violating phase (ps lH |Zl, making an im¬ 
portant contribution to the world-average value of (ps lIHHEl- The phase ps is predicted to be 
small in the standard model IITSl , making its determination interesting because of the large en¬ 
hancements that can be introduced by new physics llTSlIT^ . In what follows, we will refer to 
the fo(980) as fo and the ^(1020) as p. 

This Letter presents the measurement of the ratio of the branching fractions S(Bg —> 
]/pfo)B{fo —)■ 71+ TT^) and5(Bg —> ]/pp)B{p —> K+K^), where in both cases the J/t/^ is detected 
through its decay to The fo and the p are identified by requiring — 974MeV| < 

SOMeV and |Mk+k- “ 1020MeV| < lOMeV. The appearance of Bg —> ]/pfo decays was first 
discussed in ITT I with a theoretical estimate for Rf^/^ of approximately 0.2, which is consistent 
with results from several experiments HllllIIHUllI. Detailed studies of the tt+tt^ mass spec¬ 
trum of the Bg —> ]/p 71+ TT^ decay in 0.3 < < 2.5 GeV EUSl reveal this final state to have 

contributions from several resonances in and the fo component to range from 65.0% to 

94.5%. However, according to the same results, the contaminations from other resonances in 
IRI/t+tt- “ 974MeV| < 50MeV are several orders of magnitude lower than the fo component, 
including the non-resonant S-wave. Based on this, the measurement of Rf^/^ is performed as¬ 
suming that the selected region of is dominated by Bg —> ]/pfo decays and neglecting 

other resonances. Systematic uncertainties are assigned to the measurement owing to these as¬ 
sumptions, taking into account the uncertainty in the fo component and the interferences with 
other resonances in the selected mass window for 

Experimentally, the ratio Rig/cp is given by 

R _ 13(Bg ^ ]/pio) Bjfp -^n+TT-) _ N^^g 
io/f B{Bl^]/pp)B{p^K+K-) 

where and are the observed yields of Bg ^ with fo —> tt+tt^ and 

Bg —> ]/p{f(^ii~)p with p K+K^ decays, respectively, and efJco is the ratio of the detection 
efficiencies for the Bg decay mode with a p to the decay mode with a fo. Uncertainties in 
the b quark production cross section cancel in the ratio, as do those from the ]/p 
branching fraction and the integrated luminosity. Given the similar topologies of the two final 
states, systematic uncertainties related to the tracking efficiency and the muon identification 
also cancel in the ratio. 

2 The CMS detector 

The central feature of the GMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter. Within the 3.8 T field volume are a silicon pixel and strip tracker, a lead tungstate crystal 
electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of 
a barrel and two endcap sections. Muons are measured in the pseudorapidity range \^\ < 2.4 in 
gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid, which are 
made using three technologies: drift tubes, cathode strip chambers, and resistive-plate cham- 
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3 Event selection 


bers. Extensive forward calorimetry complements the coverage provided by the barrel and 
endcap detectors. The main subdetectors used in this analysis are the silicon tracker and the 
muon systems. 

The silicon tracker measures charged particles within the pseudorapidity range \t]\ < 2.5 and 
consists of 1440 silicon pixel and 15 148 silicon strip detector modules. Matching muons to 
tracks measured in the silicon tracker results in a relative transverse momentum resolution for 
muons with 20 < px < 100 GeV of 1.3-2.0% in the barrel and better than 6% in the endcaps. 
The px resolution in the barrel is better than 10% for muons with px up to 1 TeV ||^|. 

The first level of the CMS trigger system, composed of custom hardware processors, uses infor¬ 
mation from the calorimeters and muon detectors to select the most interesting events in a fixed 
time interval of less than 4 ps. The high-level trigger (HLT) processor farm further decreases 
the event rate to less than 1 kHz, before data storage. 

A more detailed description of the CMS detector, together with a definition of the coordinate 
system used and the relevant kinematic variables, can be found in Ref. 12T|. 


3 Event selection 

The data sample used for this measurement was collected in 2011 by the CMS experiment at 
the CERN EHC in proton-proton collisions at a centre-of-mass energy of 7 TeV and corresponds 
to an integrated luminosity of 5.3 fb^^. 

The search for Bg —>• J/t/^fo decays is performed in events with two muon candidates selected 
by the dimuon trigger at the HET, requiring the muon pair to originate from a displaced vertex. 
The dimuon candidates are further required to comply with L^y/c^y > 3, where 'Lxy is the 
magnitude of the vector Lxy, which lies in a plane transverse to the beam axis and points from 
the interaction point to the dimuon vertex, and axy is its uncertainty; cos aj/^ > 0.9, where aj/^ 
is the angle between the direction of the dimuon transverse momentum and Lxy; Pt > 4 GeV 
and Ip I < 2.2 for each muon candidate; px > 7 GeV for the dimuon; the distance of closest 
approach of each muon track with respect to the other muon track <0.5 cm. 

Reconstruction of the Bg —> J/t/^ fo decays begins with the search for J/t/? candidates by combin¬ 
ing two muons of opposite charge to form a vertex with a fit probability >0.5% and an invariant 
mass (Mj/^) within \Mj/^ — 3097.6 MeV| < 150 MeV. To search for fo candidates, two tracks of 
opposite charge assumed to be pions are constrained to a vertex with a probability >5%. One 
pion candidate must have px > 1 GeV and the other px > 2.5 GeV. In addition, the fo candidate 
must have px > 3.5 GeV and Mf^ in the range jMf^ — 974 MeV | < 50 MeV. The 974 MeV is the 
measured mass of fo signal in data modeled by a Breit-Wigner function. This value is consis¬ 
tent with the fo mass from the Particle Data Group Il22il and the EHGb measurement [IJ. Einally, 
a vertex is formed with the j/ip and fo candidates, constraining the dimuon mass to the nomi¬ 
nal ]/tp mass 11221 . The Bg ^ J/t/’fo candidates are required to have a vertex probability >10%, 
Px > 13 GeV, cos ocqo > 0.994, where ago is the angle between the direction of the Bg transverse 
momentum and the vector Ljy, and a proper decay length > 100 pm. The proper decay length 
is defined as (Lxy ■ Px Mg/Px)/ where px is the transverse momentum of the Bg candidate and 
Mg is the world-average Bg mass |22]. In the case of multiple Bg candidates per event, the one 
with smallest Bg vertex fit ^ is selected. The selection criteria for the Bg candidates are estab¬ 
lished by maximizing S / \/S -h B, where S is the signal yield obtained from Monte Garlo (MG) 
simulation and B is the background yield taken from sideband regions, defined as the number 
of events with a p^p^ invariant mass in the range 5.27 to 5.30 GeV or 5.43 to 5.46 GeV. 
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The same procedure and selection criteria are applied to the reconstruction of the normalization 
channel Bg —> ]/ip(p, except that the invariant mass requirement \M(p — 1020 MeV| < lOMeV is 
tighter than that for the fo- 


4 Results 

The signal yields of both decay channels are extracted using unbinned maximum-likelihood 
fits of the mass distributions. The invariant mass distribution of the J/t/’(^+/^^)fo(7r+7r^) can¬ 
didates is shown in Fig. It is fit with a superposition of a Gaussian function representing 
the signal, a polynomial function to account for the combinatorial background, and another 
Gaussian function for any possible peaking background. The latter models resonant structures 
that could appear in the left sideband of the signal mass owing to the 

misidentification of a kaon as a pion coming from decays such as —>• J/t/7K*(892)(K+7r“) 
and Bg —> J/t/^K+K^, as examples. In addition, B+ —> ]/xpK ^decays can be a source of 
background when combined with an extra background pion candidate. When allowing all pa¬ 
rameters to float, the fit returns = 873 ± 49 events and a Bg mass of 5369.1 ± 0.9 MeV, with 
a resolution of 15.9 ± 0.9 MeV, where the uncertainties are statistical only. The measured values 
of the Bg mass and its resolution are consistent with the MG simulation. 



J/il/TtV invariant mass [GeV] 


Figure 1: Invariant mass distribution of the Bg —)■ J/t/7(/^+/r“)fo(7r+7r^) candidates (filled cir¬ 
cles). The signal is modelled as a Gaussian (dot-dashed line), the combinatorial background as 
a first-order polynomial function (dashed-double-dotted line), and the peaking background by 
a Gaussian (dotted line). The result of the total fit is shown with the solid line. The bottom plot 
shows the pull, which is the deviation of the data from the fit divided by the uncertainty in the 
data. 

The invariant mass distribution is modelled by two Gaussian functions 

for the signal and a constant function for the combinatorial background. A signal yield of 
^obs “ 8377 ± 107 events is obtained, with a Bg mass of 5366.8 ± 0.2 MeV and a resolution 
of 17.1 ±0.1 MeV, which are consistent with the MG simulation. The corresponding invariant 
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4 Results 


mass distribution is presented in Fig. 



J/\|/K'"K' invariant mass [GeV] 


Figure 2: Invariant mass distribution of the Bg —>• candidates (black filled 

circles). The signal model is a double Gaussian (dot-dashed line), while the combinatorial 
background model is a constant function (dash-double-dotted line). The total fit is represented 
by the solid line. The bottom plot shows the deviation of the data to the fit divided by the 
statistical uncertainty in the data. 

Using the MC simulation, the detection efficiencies for the two processes are calculated as the 
ratio of the reconstructed and generated yields. The Bg meson production is simulated using 
PYTHIA 6.4.24 |[23l and its decays simulated with evtgen Il24l . The Bg mass and lifetime are 
set to 5369.6 MeV and 438/rm in the simulation. The decay model used for the Bg —> J/t/’fo 
decay is a phase-space model reweighted to reflect the spin-1 structure of the ]/rp —> 
decay. The corresponding models for the Bg —> ]/tp(p decay are: a pseudoscalar-vector-vector 
with CP violation 1251 [26l for the Bg decay, with parameters 1241 | |A|| | ^ = 0.24, | |Ao| = 0.6, 
11I p = 0.16, ^11 = 2.5, (po = 0, and (p± = —0.17; a vector-lepton-lepton model with radiation 
(photos) Il27l for the ]/xp —> decay; and a vector-scalar-scalar model 1241 for the (p —> 

K+K^ decay. The events are processed with a GEANT4-based detector simulation 1281 and the 
same reconstruction algorithms used on data. In order to validate the MC simulation samples, 
relevant kinematic and geometric variables of both simulated decay channels are compared 
with the data after background subtraction and found to be in agreement. For example. Fig. 
compares the pj and invariant mass distributions of the fo (7r+ tt” ) candidates for background- 
subtracted data and MC simulation. The fo width was set to 50 MeV in the MC simulation. 
This is consistent with what is observed in our data as shown in the Fig. |^b. The ratio of the 
detection efficiencies for the two Bg decays is calculated to be efeco = 1-344 ± 0.095, where the 
uncertainty reflects the limited size of simuated samples. Using the corresponding values of 

^obs' ^ohs' measure Rtg/cp = 0.140 ± 0.008, where the uncertainty is 

statistical only. 

The stability of the Rig/<(, measurement is verified with control checks using different run peri¬ 
ods, selection criteria, and geometric acceptances. To study possible effects from varying run 
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Figure 3: Comparison of normalized MC simulation (triangles) and background-subtracted 
data (squared) for (a) the pj and (b) invariant mass distributions of the fo (7r+ 7r“) candidates. 

conditions, the value of Rt^/ip is determined for two subsamples, found by dividing the data 
into two. The ratio is also measured after changing the selection criteria for the proper decay 
length and pj of the candidates and the pj of the leading and subleading pion candidates, 
and by using different azimuthal angle and rj requirements for the muons. None of these cross¬ 
checks revealed any statistically significant bias. 


5 Systematic uncertainties 

Potential systematic uncertainties in the measurement of Rfg/ip come from sources such as the Bg 
signal yield extraction procedure, the relative efficiency estimation, and possible contributions 
to the Bg yields from other decays producing the j/ipn^n^ and J/t/^K+K^ final states. 

Systematic uncertainties in the signal yield extraction are estimated by changing the modeling 
of the signal and the background invariant mass distributions in the likelihood fits. For the case 
of the Bg —> ]/tpio mass distribution the signal shape is changed to a double-Gaussian function 
and the background to an exponential function, while for the Bg —> ]/tp(p mass distribution 
the signal is changed to a Gaussian function and its background is modelled as a first-order 
polynomial function. These changes lead to a maximum variation of 2.1% in Rig/<f,- 

There are several factors that may affect the estimate of efeco • While the MG simulation package 
uses a Breit-Wigner model to simulate the fo —> process, it has been pointed out O |3l 

that a Flatte model is a better description of this decay. To estimate the effect of the simulation 
model, the Breit-Wigner model used in the simulation is compared to a Flatte model in the 
selected Mj^+„- region. The difference in the models reflects a systematic error of 5.8% in cfeco • 
This is quoted as a systematic uncertainty. In addition, in the MG simulation the fo width is 
set to 50MeV. This value is varied by ilOMeV, resulting in a systematic uncertainty of 8.6% 
in Rfp/^. The models used in the MG simulation of the Bg decays are set to phase-space 12411 
instead of the default decay models, leading to a 6.2% systematic uncertainty in Rfg/^. Finally, 

the statistical uncertainty in cfeco owing to the finite number of MG events, which corresponds 
to 7.1%, is added as a systematic uncertainty. 

As mentioned in the introduction, detailed studies of the tz^ tz^ mass spectrum of the Bg —> 
]/\p TZ^TZ^ decay 121131 in a mass window of 0.3-2.5 GeV, reveal this final state to have contribu- 
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6 Summary 


tions from several resonances in and the fo component to range from 65.0 to 94.5% in 

the entire mass window studied by LHCb. To study the effects of the interferences and the fo 
fraction observed by LHCb in the estimate of efeco / the model reported in ||3l for the lowest fo 
fraction and largest non-resonant component was compared to the single Breit-Wigner model 
used in the MC simulation of the fo —> tt+tt decay. The comparison in the selected re¬ 

gion shows a variation of 5.6% in This is quoted as a systematic uncertainty coming from 
this source. It can be observed in the same LHCb study that the contaminations from other 
resonances in the mass region — 974MeV| < 50MeV are several orders of magnitude 

lower than the fo component, including the non-resonant S-wave. To estimate the variation in 
the Bg yield coming from these possible contributions, the fo mass window is widened from 
50 to 100 MeV around the fo mass, resulting in a variation in Rf^/^ of 6.4% that is quoted as a 
systematic uncertainty. For the Bg —)■ J/i/?K+K^ decay charmel, the contribution of the S-wave 
in a ^ mass window similar to what is used in this analysis has been found to be negligible Il29l . 

Combining these uncertainties in quadrature leads to a total systematic uncertainty of 16.5%. 


6 Summary 

Using data collected by the CMS experiment in proton-proton collisions at y/s = 7TeV, corres¬ 
ponding to an integrated luminosity of 5.3fb^^, 873 ± 49 events of Bg —> J/t/?(f/+f/“)fo(7r+7r“) 
and 8377 ± 107 events of Bg —>• J/i/?(/r+^^)^(K+K^) are observed. The fo and (p are identi¬ 
fied in the mass ranges — 974MeV| < 50MeV and |Mk+k- “ 1020MeV| < 10MeV, 

respectively. The ratio of the branching fraction of Bg —> J/t/;(^+f/~)fo(7r^7r^) to the branching 
fraction of Bg —> is found to be 


g(B° ^ g(fo ^ 71+71 ) 

6 (BO ^ ]/tp(p) B{(p^ K+K-) 


0.140 ± 0.008 (stat) ± 0.023 (syst). 


( 2 ) 


This result is consistent with the theoretical prediction of about 0.2 IITtI and with previous 
measurements in different ranges of EIHIl^. 
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